Analytical study of quadratic and non-quadratic 
short-time behavior of quantum decay 
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The short-time behavior of quantum decay of an unstable state initially located within an inter- 
action region of finite range is investigated using a resonant expansion of the survival amplitude. It 
is shown that in general the short-time behavior of the survival probability S(t) has a dependence 
on the initial state and may behave either as S(t) = 1 — 0(t ' ) or as S(t) = 1 - 0(t 2 ). The above 
cases are illustrated by solvable models. The experiment reported in Ref. [1] does not distinguish 
between the above short-time behaviors. 
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Introduction. The decay of unstable systems, corre- 
sponding to particle emission by tunneling out of a poten- 
tial, has been a subject of attention since the early days 
of quantum mechanics. In 1928 Gamow derived an ex- 
pression for the exponential decay law and introduced the 
notion of the lifetime [2], which provides the time scale 
for exponential decay and sets the meaning for short or 
long times in decay. In the fifties of last century, Khalfin 
demonstrated that if the energy spectra of the system is 
bounded by below, the exponential decay law cannot hold 
at long times [3]. At short times there is also a depar- 
ture from the exponential decay behavior which, however, 
it is related to the energy moments of the Hamiltonian 
to the system [4-7]. It is usually assumed that decay 
at short-times exhibits a quadratic behavior with time 
[4, 5, 7-10]. The relevant quantity is the survival prob- 
ability S(t) = \A(t)\ 2 , with A(t) = (0|exp(-iHf/ft)|0>, 
that yields the probability that at time t the system re- 
mains in the normalized initial state |0). Notice that 
expanding exp(— itlt/h) yields 
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which leads to an expansion of S(t) that involves only 
even powers of t. A quadratic behavior requires that at 
least the first two moments of the Hamiltonian to the 
system are finite. The experimental verification of the 
short-time behavior of decay was provided some years 
ago and seems to be consistent with an initial quadratic 
behavior [1]. However, from a theoretical point of view, it 
is not obvious that the series expansion of S(t) mentioned 
above converges or even that it is defined. In this context, 
there are problems that have been rarely explored as the 
conditions that may lead to a non-quadratic behavior at 
short times [7, 11]. Some recent work has discussed a i 3 / 2 



short-time behavior of decay in the context of specific 
models [12-15]. It is also worth mentioning that a t 3 / 2 
short-time behavior has been found in studies involving 
transients in non-decay problems [16]. 

Here we consider an approach to the time evolution 
of decay based on a resonant expansion of the survival 
amplitude that has been studied intensively for the expo- 
nential and long-time regimes [17-19]. The occurrence, 
however, of a double sum in the expression for S(t), which 
in general does not commute, prevented its application 
to the discussion of the short-time behavior of decay [20] . 
Recently, however, motivated by the considerations men- 
tioned in the previous paragraph, we believe that we have 
found a way to circumvent the above situation. 

In this work we address a rigorous investigation of the 
short-time behavior of decay for unstable systems. We 
obtain general conditions on the initial states so that S(t) 
may exhibit either a time dependence as t 3 / 2 or as t 2 . 
We also indicate that the experiment on the short-time 
decay behavior reported in Ref. [1] does not distinguish 
between the above two cases. 

Resonant expansion. We consider a simple yet no triv- 
ial description of the decay process that involves real po- 
tentials of arbitrary shape that vanish beyond the interval 
{0,L), which is well justified since most effective poten- 
tials in physics are of short-range, and initial states that 
are confined initially within the interaction region. The 
above conditions are commonly found in quantum sys- 
tems designed artificially, as low temperature multibar- 
ricr resonant tunneling structures [21] or ultracold atoms 
confined in optical traps [22]. A relevant feature of these 
systems is that the decay process is essentially coherent 
(clastic). One may then exploit the analytical properties 
of the outgoing Green's function to the problem on the 
whole complex wave number plane where it possesses an 
infinite number of poles. These poles are in general sim- 
ple and are distributed in a well known manner [23] . This 
has led to a formulation of the time evolution of decay 
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in terms of a purely discrete expansion that involves the 
residues (resonant states) at the poles of the outgoing 
Green's function to the problem [17-19]. The resonant 
states u n (x) satisfy the Schrodingcr equation of the prob- 
lem obeying purely outgoing boundary conditions and 
hence they also include the bound and antibound states 
of the problem. 

It is worth stressing that the resonant state formula- 
tion yields exactly the same results as a calculation using 
continuum states [24]. 

One should mention that the above analytical proper- 
ties of the the outgoing Green's function remain valid for 
potentials having tails that go faster than an exponential 
at infinity, as for example having Gaussian tails [23] . The 
outgoing Green's function for potentials having exponen- 
tial tails may be extended analytically only through a 
finite region of the complex k plane, and hence the corre- 
sponding expansion will consist in addition to a discrete 
pole expansion of an integral contribution involving con- 
tinuum of states. We believe, however, that this issue 
is mostly of mathematical interest since as pointed out 
above, most effective potentials in physics are negligible 
small after a distance and hence are beyond experimental 
scrutiny. 

The survival amplitude may be expanded in terms of 
resonant states as [17-19] 
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where u(iy n ) refers to the Faddeyeva function [25] with 
y n = - exp(-in /4)(h/2m) 1/2 Knt 1 / 2 , and n n = a n - i/3 n 
relates to the complex energy eigenvalue E n = h 2 K 2 l /2m. 
Notice that for bound states n n = i"/ n with 7„ > 0, 
and similarly for antibound states with 7„ < 0. The 
coefficients C n and C n in (2) are 

C n = / ip (x, 0) u n (x) dx; C n = / ip* (x, 0) u n {x) dx. 
Jo JO 

(3) 

The above coefficients fulfill the relationship [18, 19] 
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Notice that Eq. (2) follows using 

OO 

yj{x,t) = (1/2) C n u n (x)cj(iy n ), (7) 



where < x < L [17-19]. Since w(0) = 1 [25], then 
^(x,0) = (l/2)X)£L-oo C n u n [x). Using that Hu n (x) 
/.',:'■'• ' and the definition of C n given in (3), allows to 
express the first moment of H as 
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which is a finite quantity as follows by inspection of the 
conditions satisfied by the potential and the initial wave 
function. 

The function u>(iy n ), which may be evaluated by well 
developed numerical methods [26] , may be written as the 
convergent expansion (for any value of iy n ) [25], 

^) = E r(it/2) ' a = e-^(h/2mf 2 . (9) 

One may write, therefore, uj(iy n ) at short times as 
w(iy n ) = 1 + aK n t x l 2 + a 2 K^t + . . .. Substitution of this 
expression into (2) allows to write A(t) at short times as 
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The term with s = 3 in Eq. (10) reads 

3 / 00 \ 

C n C n K n )t 3 / 2 . (11) 
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However, depending on the characteristics of the initial 
state ?p(x,0), the sum in (11) may vanish, be a constant 
or diverge. Let us first analyze the case where it vanishes. 
Then, we may write (10) with the sum over s starting 
from ,s = 4 in the form 
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A(i) = l--(H)i+- ]T C n C n O n (t 2 ), (12) 



where we have used, respectively, Eqs. (4), (5), (8) and 
O n (t 2 ), with O the O-symbol [27], expresses the fact that 
as t — > 0, the leading term in the remaining absolutely 
convergent sum over s is t 2 . Hence the survival probabil- 
ity may be written as 



s(t) = i + -(u) 2 t 2 -o(n, 



(13) 



Since as t ->• 0, O n (t 2 )lt v 0, requires that v < 2, it fol- 
lows that 0(V)/t u — > provided r > 2. Notice however 
that since 5(0) = 1, the decay process implies S(t) < 1 
for t > 0. Hence, in order to avoid that the term pro- 
portional to t 2 in (13) yields an unphysical time interval 



3 



where S(t) > 1, necessarily r — 2. This guarantees that 
S(t) diminishes with time and hence 



S(t) = 1 - Oif 



(14) 



The above result seems to hold independently of whether 
or not the second moment (H 2 ) oc Y^=-oo CnCnK^ is fi- 
nite; the second case, where the sum in (11) is a constant, 
gives that the leading term of S(t) at short times is t 3 ^ 2 , 
and the last case, where the sum in (11) diverges, implies 
that such a term cannot be extracted from the sum over 
s in (10), i.e., for each value of n one has to perform the 
convergent sum over s, and hence 
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Notice that otherwise it would be a contradiction with 
the argument leading to Eq. (14), that rests on the as- 
sumption that Eq. (11) vanishes exactly. Hence neces- 
sarily r = 3/2, and therefore as t — > 0, 



S(t) = 1 - 0(t 3/2 ) 



(16) 



It is worth mentioning here that in Ref. [7] reports the 
possibility of a t 3 ! 2 short-time dependence of S(t) pro- 
vided the second energy moment diverges. Though ref- 
erence to a pole expansion of A(t) involving w-functions 
is made to account for a possible fractional behavior, the 
analysis there is actually based on the finiteness or not 
of the expressions 
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[A(t)] t=0 = -(i/fi)(H) 



[A(t)] t=0 = -(l/h 2 )(H 2 ), 
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the dot indicating derivative with respect to time, which 
are obtained from the scries expansion of exp(— iHt/h). 
Using the exact expansion (2) one sees that [A(t)] t= o re- 
mains as given by (17) above, but 

oo ^3 oo j 

A(t) = a 4 2J C n C n n 4 i uj(iy n ) + — = ^ C n C n i£ n -^ 
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(19) 

is different. One sees therefore that [-A(t)]t=o diverges 
unless (11) vanishes, independently of whether (H 2 ) is 
finite or not, contrary to the result given in Ref. [7]. 
This indicates that the series expansion of exp(— iHt/H) 
may lead to misleading results. 

The short-time expressions given by Eqs. (14) and (16) 
suggest to consider the expression for S(t) 



S(t) 
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with parameters i? and t* , to adjust the short-time be- 
havior of exact calculations, using (2), or experiment. 



Model. As pointed out above, the distinction between 
the t 3 / 2 and t 2 short-time behavior of S(t) depends on the 
properties of the initial states. Theoretically, this nec- 
essarily leads to model calculations. For simplicity, we 
consider as a model a double-barrier resonant tunneling 
nanostructure [21, 28] that extends from to L = 15 nm, 
with 6 = 5 nm (barrier widths), w = 5 nm (well widths), 
V = 0.23 eV (barrier heights) and m = 0.067m e (effec- 
tive electron mass). There exist well developed proce- 
dures to obtain the set of resonant states {u n (x)} and 
complex poles {n n } for a given problem [19, 29]. We 
choose two different types of initial states within the in- 
teraction potential region < x < L: a cutoff Gaussian 
pulse 



ip(x) = (l/27r t 7 2 ) 1/4 exp[-(x-.T ) 2 /4cr 2 



(21) 



centered at x$ = L/2 with pulse width a <C w, to guar- 
antee that the effect of cutting-off the tails is negligible, 
and a sinusoidal pulse: 



ipj(x) = \j2jw sin [kj (x — b)] 



(22) 



with b < x < b + w, and zero elsewhere, where kj = jir/w 
for a fixed integer value j. 

The reason for the above choice of initial states is 
that, in addition to mathematical simplicity, one ex- 
pects on physical grounds that the decaying particle 
is initially confined within the interaction region and 
hence that possible tails beyond that region are negligi- 
ble. Of course, one may envisage an initial state having 
large non-negligible external tails. In that case, as time 
evolves, part of the external portions of the initial state 
would head towards the internal region and would inter- 
fere with the decaying part giving origin to a transient 
behavior. We are not addressing such a possibility in this 
work though it might be of interest to investigate it. 

In order to study numerically the behavior of the sur- 
vival probability at short times, it is convenient to define 
the quantity 



A N (t) 
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A' 
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Hence, the moments of the Hamiltonian may be writ- 
ten as (W) = (l/2)(h 2 /2m)nim N ^ 00 n N (2j). Since 
A*aK0) — > 2 when N — > oo, one finds that 5jv(i) = 
(20) \A N (t)\ 2 fulfills lim^oo S N (t) -> S(t). Here, S N (0) = 1 
for any value N > 0. Moreover, if for two values N 
and N' such that N N' the corresponding survival 
probability satisfies Sjv(i) = S^'it) in a time interval, 
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FIG. 1. (color online) Plot of In Sat (t) as function time in 
lifetime units at short times for (a) An initial Gaussian pulse 
and (b) An initial sinusoidal pulse. Both cases exhibit a com- 
parison of the exact calculation (2), respectively, for N = 10 3 
(dotted line) and N = 2 x 10 4 poles (solid line). Each figure 
exhibits also the corresponding comparison with (20) (dashed 
line), with adjustment parameters ■& — 2 and r* = 0.819 fs in 
(a), and i9 = 3/2 and r* = 3.802 fs in (b). See text. 



that implies that both approximations yield the correct 
behavior of S(t). 

We have evaluated Eq. (11) for N = 10 3 and N = 
2 x 10 4 for both the Gaussian and the sinusoidal pulses 
and found that it vanishes for the Gaussian pulse and di- 
verge for the sinusoidal pulse. This implies, according to 
our analysis above, that the initial Gaussian state should 
produce a quadratic short-time behavior whereas the ini- 
tial sinusoidal state a non-quadratic one. This needs to 
be confirmed by a comparison between an exact calcu- 
lation, using (2), and the adjustment formula (20). We 
have also evaluated /ijv(j) for values j = 4, ...,8 for the 
above pulses and found that these quantities are finite for 
the Gaussian pulse and diverge for the sinusoidal one. 

Figures la and lb, exhibit, respectively, plots of 
lnSW(t) at short times in units of the lifetime t%, for the 
initial Gaussian pulse with a = w/10 and for the initial 
sinusoidal pulse with j = 1, using the exact pole expan- 



sion (2) with N = 10 3 poles (dotted line) and N = 2x 10 4 
(solid line). One sees, in each figure, that these curves are 
indistinguishable from each other which indicates excel- 
lent convergence using N — 10 3 poles. Each of the above 
figures also exhibits the results of the calculation employ- 
ing the adjust formula (20) employing, respectively, the 
origin and two other points of the corresponding exact 
calculation. In Fig. la the adjustment yields $ = 2 and 
t* w 0.819 fs (dashed line) which confirms the quadratic 
short-time behavior given by Eq. (14). In this case, we 
find that r* « tz where tz is the Zcno time defined by 
t z = h/AE, with AE 2 = (H 2 ) - (H) 2 [10]. We have ob- 
tained similar results for initial Gaussian states having 
different values of the width a provided a « ic. Sim- 
ilarly, in Fig. lb, the adjustment (dashed line) yields 
i9 = 3/2 and r* = 3.802 fs which confirms the fractional 
t 3 / 2 short-time behavior given by Eq. (16). Similar re- 
sults occur for other values of j. 

It is worth emphasizing that if the first two moments 
of H exist in the expansion given by Eq. (1), then con- 
sistency with the expansion given by Eq. (2), which in 
general involves both quadratic and non-quadratic pow- 
ers of t, requires that the term proportional to t 3 / 2 should 
vanish, as indeed we have numerically corroborated. 

Experiment. The experiment of Rcf. [1] involves an 
external potential that goes linearly with distance and 
hence our analysis is not strictly applicable. However, 
we find of interest to perform an elementary adjustment 
using (20) to the data given in Rcf. [1], which assumed a 
quadratic short-time behavior [1, 30]. Since our analysis 
predicts that the value of $ in Eq. (20) is either 2 or 
3/2, we need to consider only two experimental points 
to make the adjustment. We choose the points with a 
minimum error bar, in particular at t = 0, and these 
correspond to Fig. 3b of Ref. [1]. In Fig. 2, we plot the 
experimental data of Fig. 3b at very short-times and the 
corresponding adjustments using Eq. (20) for $ = 2 and 
t* = 12.55 /is (full line) and = 1.5 and r* = 23.15 
(dotted line). We see that both short-time behaviors are 
consistent with experiment yet with a different value of 
the time scale r*. May be a future experiment could 
discriminate between these two time scales. 

Concluding remarks. It is worth emphasizing that in 
general, the expansion of A(t) in powers of t is not de- 
fined. This means that the corresponding Taylor expan- 
sion around t = does not exist in general. The van- 
ishing or not of Eq. (11), which determines a quadratic 
or non-quadratic time evolution at short times, is very 
sensitive to the tails of the initial state, as exemplified 
by the Gaussian and sinusoidal initial states discussed 
here. It is not clear, therefore, that initial physical states 
possess finite moments, a point that has been a subject 
of debate [31]. Further study on the characterization of 
initial states is needed [32]. It is also worth mention- 
ing that a non-quadratic £ 3 / 2 short-time behavior docs 
not prevent the occurrence of the quantum Zeno effect 
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FIG. 2. (color online) Plot of the numerical adjustments at 
short times, using Eq. (20), to the corresponding experimen- 
tal data taken from figure 3b of ref. [1]. See text. 

[5, 13, 15]. Finally, we believe that our results may be 
relevant for quests regarding the description of the short- 
time behavior of unstable systems in relativistic quantum 
field theory where it has been found that the second mo- 
ment to the Hamiltonian diverges [34]. 
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